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Abstract 

The vacuum polarisation effects from tlie nucleon sector lead to large medium modifications of 
the vector meson masses in the Walecka model. With a quantum hadrodynamic framework includ- 
ing the quantum effects, and using a quasiparticle description for quark gluon plasma(QGP), the 
dilepton emission rate from the hot and dense matter resulting from relativistic nuclear collisions 
is calculated. Using a model for the fireball evolution which has been shown to reproduce other ob- 
servables such as charmonium suppression, photon emission and the abundance of hadronic species, 
we compare with the dilepton invariant mass spectrum measured by the CERES collaboration at 
CERN SPS. We also compare the results to a previous calculation where the spectral function of a 
virtual photon (a key ingredient for the emission rate) has been calculated in a chiral SU(3) model. 
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1. INTRODUCTION 



The study of the in-medium properties of the vector mesons (p and u) in hot and dense 
matter is actively investigated, both experimentally [l[ 0| and theoretically 0, Q]. The 
experimental observation of enhanced dilepton production [l^ in the low invariant mass 
regime possibly is due to a reduction in the vector meson masses in the medium. It was first 
suggested by Brown and Rho that the vector meson masses drop in the medium according 
to a simple scaling law jsf. Within the framework of Quantum Hadrodynamics (QHD), the 
vacuum polarisation effects from the baryon sector lead to a significant drop of 

the vector meson masses in the medium whereas the mass modification is marginal with only 
Fermi sea polarisation effects. The properties of the hadrons as modified in the thermal bath 



are reflected in the dilepton and photon spectra emitted from a hot and dense matter 
0, ^|. Dileptons are interesting probes for the study of the hot and dense matter formed 
in relativistic heavy ion collisions. Since they do not interact strongly, the dileptons escape 
unthermalized from the hot and dense matter at all stages of the evolution. The observed 
enhancement of dileptons in the low invariant mass regime has initiated extensive theoretical 
investigations on the temperature and density [1^ modifications of the dileptons from 
hot hadronic matter as well as from a quark gluon plasma (QGP) formed in heavy ion 
collisions. The increased dilepton yield in S+Au collisions as observed by the CERES 
collaboration was attributed to enhanced p-meson production via tt^tt^ annihilation and a 



dropping of the p-mass in the medium \l4 . Il5l | . A lar ge b roadening of the p-meson spectral 
function arising due to scattering off by baryons [la ll2| has been shown to reproduce the 
CERES data quite well. In the present investigation, the space-time evolution of the strongly 
compressed hadronic matter formed in a relativistic heavy ion collision is considered in a 
mixed scenario of QGP and hadronic matter. The QGP is described within a quasiparticle 
picture and the hadronic matter described in a Quantum Hadrodynamic framework 
including quantum fluctuation effects from the baryon and scalar meson sectors 
The dilepton emission rates from the expanding fireball are then studied in a model based 
on local thermal equilibrium and isentropic expansion 2^. 

It was earlier demonstrated in a nonperturbative formalism that a realignment of the 
ground state with baryon-antibaryon condensates is equivalent to the relativistic Hartree 
approximation (RHA) |20|. The ground state for the nuclear matter was extended to include 
sigma condensates to take into account the quantum correction effects from the scalar meson 
sector 0. Such a fonnal.m includes mu.tUoop effects aud is self cous,stent HQ. The 
methodology was then generalized to consider hot nuclear matter [2^ as well as to study 
hyperon-rich dense matter j^l relevant for neutron stars. The effect of vacuum polarisations 
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22 1 has also been 



on the vector meson properties j21[ and on the static dilepton spectra 
recently studied. In the low invariant mass regime, the scalar meson contributions lead to 
considerable broadening the u peak in the dilepton spectra as compared to RHA, which leads 
to smearing and ultimate disappearance of the uj peak at high densities 2^. The present 
investigation of dilepton spectra using a dynamical fireball model is compared with results 
obtained using other descriptions for hadronic matter jl^ . \v\ , as well as to the experimental 
results from CERES collaboration 

We organize the paper as follows. We first briefly recapitulate the nonperturbative frame- 
work used for studying the medium modification of the uo and p vector meson properties 
including the vacuum fluctuation effects in the Walecka model, in section 2. Section 3 dis- 
cusses the parametrisation of the photon self energy in the hadronic phase. The quark gluon 
plasma (QGP) phase is treated within a quasiparticle picture [3] as described in section 
4. In section 5, we describe briefly the model considered for the dynamical evolution of the 
fireball. Section 6 deals with the calculation of dilepton emission from the fireball. Finally, 
we discuss the findings of the present investigation in section 7 and summarise in section 8. 



2. VACUUM POLARISATION EFFECTS AND IN-MEDIUM VECTOR MESON 
PROPERTIES 

We briefly recapitulate here the vacuum polarisation effects arising from the nucleon and 
scalar meson fields in hot nuclear matter in a nonperturbative variational framework [2^ 
and their influence on the vector meson properties in the hot and dense matter. The method 
of thermofield dynamics (TFD) is used here to study the "ground state" (the state with 
minimum thermodynamic potential) at finite temperature and density within the Walecka 
model with a quartic scalar self interaction. The temperature and density dependent baryon 
and sigma masses are also calculated in a self-consistent manner in the formalism. The ansatz 
functions involved in such an approach are determined through functional minimisation of 
the thermodynamic potential. 

The Lagrangian density in the Walecka model is given as 

+ \mluj^'uj, - -^{d^u, - d,uj^){d^'u'' - d^u^). (1) 

In the above, if), a, and uj^ are the fields for the nucleon, ex, and uj mesons with masses M, mg-, 
and respectively. The quartic coupling term in a is necessary for the sigma condensates to 
exist, through a vacuum realignment [2^. Our calculations thus include the quantum effects 
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arising from the sigma meson in addition to the mean field contribution from the the quartic 
self interaction of the scalar meson. We retain the quantum nature of both the nucleon and 
the scalar meson fields, whereas the vector uj- meson is treated as a classical field, which in 
the mean field approximation is given as (cj^) = (5^o^o- The reason is that without higher- 
order term for the cu-meson, the quantum effects generated due to the cu-meson through the 
present variational ansatz turn out to be zero. The thermodynamic quantities including 
the quantum effects can then be written down. The details regarding the formalism can be 



found in earlier references 



. l25| . Extremisation of the thermodynamic potential, with 
respect to the meson fields cxo and ojq yields the self-consistency conditions for do (and hence 
for the effective nucleon mass, M* = M + g^ao) , and for the vector meson field c^o- We might 
note here that the quantum effects arising from the scalar meson sector through a meson 
condensates amount to a sum over a class of multiloop diagrams and, do not correspond to 
the one meson loop approximation for scalar meson quantum effects considered earlier . 



2.1. In- medium vector meson masses 

We now examine the medium modification to the masses of the u- and p-mesons in hot 
nuclear matter including the quantum correction effects in the relativistic random phase 
approximation. The interaction vertices for these mesons with nucleons are given as 



Cr^t = 9v (^7m^>K^ - (2) 

where VJ^ = uo^ or p^, is the free nucleon mass, is the nucleon field and = 1 or r, 
r being the Pauli matrices, gy and ny correspond to the couplings due to the vector and 
tensor interactions for the corresponding vector mesons to the nucleon fields. The vector 
meson self energy is expressed in terms of the nucleon propagator, G{k) modified by the 
quantum effects. This is given as 

W^{k) = -li9vj^, I d'pTT[T>i.{k)G{p)TU-k)G{p + k)\, (3) 

where 7/ = 2 is the isospin degeneracy factor for nuclear matter, and Ty^k) = I'^'^a — ^^cr^" 
represents the meson-nucleon vertex function. For the u meson, the tensor coupling , being 
small as compared to the vector coupling to the nucleons , is neglected |2l[ I2J. After 
carrying out the renormalization procedures for the vector self energies, the effective mass 
of the vector meson is obtained by solving the equation 



kl-ml + Ren(A;o, fc = 0) = 0. (4) 



where 11 = ^Uf,. 
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2.2. Meson decay properties 



For a baryon-free environment, p — > tt+vt" is the dominant decay channel for p meson. 
The decay width for this process is calculated from the imaginary part of the self energy 
using the Cutkosky rule, and in the rest frame of the p-meson is given by 



487r 



^2 



(5) 



where, f{x) = [e^^ — is the Bose-Einstein distribution function. The first and the second 
terms in the equation correspond to the decay and the formation of the resonance, p. In 
the calculation for the p decay width, the pion has been treated as free, and any modification 
of the pion propagator due to effects like delta-nucleon hole excitation to yield a finite 
decay width for the pion, have not been taken into account. The coupling gp^^T^ is fixed from 
the decay width of p meson in vacuum (rp=151 MeV) decaying to two pions. 

In the presence of baryons, however, it has been shown that there is considerable increase 
of the p decay width in the thermal medium due to the scattering off the baryons. The 
dominant contributions which lead to appreciable broadening of the p- spectral function are 
the inelastic processes pN —>■ ttN and pN — > vrA. The imaginary parts of the corresponding 
scattering amplitudes in the large baryon mass limit are given as jlT^ 



ImT 



(ttN) 
pN 



g\n{ko,0,m^) 



and 
where, 



ImT 



(^A) 
pN 



'^gAH{ko,m*^ 



(6) 



(7) 



gpNN 



^(ko-Ay-m^ 



{3k^ - Aklm^ + 4m^ + A(8fcom^ - I2kl) 



+ A^{lQkl-Sm^) 



8A^ko + 4A^ 



{3kl 



4A^ 



(8) 



2A;o(A;2-4A2)2 

In the above, gA_ is the axial vector coupling constant chosen to be gA = 1-26 17J. The pNN 
coupling is as fitted from the A^A^ scattering data 30] and is given as gpNN='^-Q- We might 
note that, by definition, the pNN coupling used here is half of the coupling in Ref |12(. 

For high energies, the pN scattering process, pN —>■ uN also becomes important jivj l. 
The corresponding scattering amplitude is given as 



ImT, 



ujpn 



2k^ Tn*j^ 



mlf{k^o - ko)[{k^o - - ] 



pN 



487r/2m2 m^(m^ + /i;o)(m2 



m* 2 + 2k^oko 



klY 



(9) 
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The interaction Lagrangian describing tupvr as used above, is given as [32| 

Hu^p. = ^e.^apd'^Lu'dy^TT, (10) 

For the iopir couphng we take the value (yf^p^ =2 according to Ref. jsl consistent with 
the decay width of ^ 717. The ujpir couphng as here is similar to that used in 
after accounting for a factor of m,r//7r difference in the definitions. In the above, {k^^l = 

A(m^ + A;o, , m^) / (2(m^ + /eq) and ^^^0 = \/ m* ^ + /c^^ . The Kallen function A is 
defined by \{x^y,z) = (x^ — + zy){x'^ — iu — 2)^)- The present calculations for the 
collisional decay width of p meson take into account the mass modifications of the nucleon 
and vector mesons (p and u). For A, we assume that the mass scales in the same way as the 
nucleon mass in the medium, with mA=1232 MeV as the mass in vacuum. It may be noted 
that the contributions given by (jB)), ((Zj) and to the pN scattering amplitude correspond 
to the zero temperature situation. The broadening of the p spectral function is significantly 
due to the finite density effects, which are expected to dominate over the temperature effects 
for the SPS conditions 01 • 

The contribution to the decay width of p due to pN scattering processes, arising from the 
ttN, ttA and ujtt loops are given in terms of the imaginary parts of the scattering amplitudes 
as 

coll yBK-^'X'-^pN -r^""p7V ^^'"'-'pN 



pN = z y^^) 



ko 

We might note here that the effect of pN scattering on the real part of the scattering 
amplitude has been seen to be negligible and has not been considered here. 

To calculate the decay width for the o^-meson, we write down the effective Lagrangian 
for the u meson as 
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2 

= -^u^A, + ^e^,^^e,,kUJ^d''7r'd''7r^df7r\ (12) 

In the above, the first term refers to the direct coupling of the vector meson uj to the photon, 
and hence to the dilepton pairs, as given by the vector dominance model. The decay width 
of the cu- meson in vacuum is dominated by the channel u —>■ Sn. In the medium, the decay 
width for Svr is given as 

= ^ j d%d%d%5^''\P-pi-p2-p3)\Mf,\^ 

;i + /(Ei))(l + /(E2))(l + f{Es)) - f{E,)f{E2)f{Es)] , (13) 

where d^pi = j2^i^, Pi and Ei's are 4-momenta and energies for the pions, and /(i?j)'s are 
their thermal distributions. The matrix element Mfi has contributions from the channels 
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pn Sn (described by Eq. (fTUj) and the direct decay u ^ Sir resulting from the 



contact interaction (second term in (fT^ l 



32, 



331 . 134| . After fitting g^p.„ from the u ^ tcj 
decay width, the point interaction couphng (7^3^ is determined by fitting the partial decay 
width C(j ^ Stt in vacuum (7.49 MeV) to be 0.24 22|. The contribution arising from the 



direct decay is up to around 15 %, which is comparable to the results of j32, |3 

With the modifications of the vector meson masses in the hot and dense medium, the 
contribution fro,„ the decay mode . ^ becomes accessible B0 when < > + 
This is also considered in the present work. There may also be collisional broadening effects 
jl7l | for the u meson at high eriergies due to baryons. However, these are seen sensitive to 
the meson-baryon form factors [l7[ and have not been considered in the present calculations. 



3. PARAMETRISATION OF THE PHOTON SELF ENERGY IN THE 
HADRONIC PHASE 

Lattice simulations indicate that QCD undergoes a phase transition at a critical temper- 
ature, Tc, above which there is QGP phase, with quarks and gluons as the relevant degrees 
of freedom. As one approaches the critical temperature from the above, the quarks and 
gluons get confined to form the hadrons, which become the effective degrees of freedom for 
temperatures below Tc. In this section, we discuss the parametrisation of the photon self 
energy in the present scenario for the hadronic phase. The photon here couples to the vector 
mesons (p and 00) and the electromagnetic current-current correlator can be related to the 
currents generated by these mesons which are calculated using a Lagrangian describing the 
hadronic phase. The averaged photon spectral function is defined as 0] 

1 2tt 

R{q) = — Imn(g), (14) 
where the photon self energy is related to the sum of the vector meson self energies as 

lmniq) = -\j2nf{q). (15) 

^ V 

In the limit of 0, the averaged photon spectral function can be written as 

R{q) = 127rg2^Imn^(go) = 127r ^ Imn^(go). (16) 

V V 

In the Walecka model, the vector self energies can be parametrised as a Breit-Wigner 
distribution with an energy dependent decay width, along with a continuum as 

^2 qo^viqo) , Cvf^ , as(go)^ 1 



7 



for the vector mesons, V = p,uj. In the above, fp and are the couphng constants 
corresponding to the vector- meson-photon interactions. The continuum part is described 
as a smooth function in the en ergy IsgI IstI i instead of a step function as is usually adopted 



in QCD sumrule calculations [17|, |38|. The values, Cp=l and C^=l/9 give the correct 



asymp totic limits for the spectral functions, limg^^co Rp — 3/2 and limgg_^oo =1/6 [1 



37 



38| . Also, the energy dependence of the strong coupling constant occurring in the meson 
self energies has been taken into consideration and is given as as{qo) = 0.7/ ln(go/0.2) jsfij ]. 
where go is in units of GeV. The values for the parameters q^'^ correspond to the threshold 
energies above which asymptotic freedom is restored and quark model estimates for cross- 
sections become valid. The vacuum values for the parameters appearing in the continuum 
part of the spectral function as derived from the experimental data of e^e~ hadrons 
are = 1.3 GeV, q^ = 1.1 GeV, and 5 = 0.2 GeV The vacuum values for the other 
parameters in the Breit-Wigner part of the photon spectral function are /p=152 MeV, /a;=50 
MeV, mp=770 MeV, m<^=783 MeV, 1^=151 MeV, =7.5 MeV. The masses of the vector 
mesons are replaced by the medium modified masses including quantum correction effects 
|22| |. as a simple pole approximation However, the energy dependence of the decay 

widths is considered. We assume the medium modified threshold energies, q^* to be given 
by the simple scaling [s^ = The above parametrisation of the photon spectral 

function for the hadronic phase is used for studying the dynamical evolution of the strongly 
interac ting matter arising in a ultrarelativistic heavy ion collision in an expanding fireball 
model [2J|. In the following section, we briefly discuss the parametrisation of the photon 
self energy in the QGP phase described by a quasiparticle picture jioj ]. 



4. PHOTON SELF ENERGY IN THE QGP PHASE 

The QGP phase is described [3| using a quasiparticle picture. The model treats quarks 
and gluons as massive thermal quasiparticles with their properties determined so as to 
be compatible with the lattice QCD data. For temperatures much larger than T^, the 
thermal masses can be calculated in the hard thermal loop (HTL) approximation using 
the perturbative QCD techniques. For temperatures near the phase transition, however, 
the coupling as becomes large thus invalidating any perturbative techniques. Around Tc, a 
power law fall-off is assumed for the thermal masses, based on the conjecture that the phase 
transition is either weakly first order or second order as indicated by the lattice calculations. 
The thermodynamic quantities for the QGP are calculated in terms of two functions B{T) 
and C{T) introduced in the model, which account for the thermal vacuum energy and the 
onset of confinement for temperature approaching Tc. The quasiparticles are by construction 
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noninteracting and hence the quasiparticle qq becomes the only contribution to the photon 
self energy |l9|. In the QGP phase the time like photon couples to the continuum of thermally 
excited qq states and subsequently converts into a dilepton pair, which gives the contribution 
to the dilepton emission rate from the QGP phase. 



5. THE FIREBALL EVOLUTION MODEL 
5.1. Expansion and flow 

Our fundamental assumption is to treat the fireball matter as thermalized from an initial 
proper time scale tq until breakup time tj. For simplicity, we assume a spatially homogeneous 
distribution of matter. Since some volume elements move with relativistic velocities, it is 
sensible to choose volumes corresponding to a given proper time r for the calculation of 
thermodynamics, hence the thermodynamic parameters temperature T, entropy density s, 
pressure p, chemical potentials /ij and energy density e become functions of r only for such 
a system. In the following, we refer to r as the time measured in a frame co-moving with a 
given volume element. 

In order to make use of the information coming from lattice QCD calculations, we proceed 
by calculating the thermodynamical response to a volume expansion that is parametrized in 
such a way as to reproduce the experimental information about the flow pattern and HBT 
correlations as closely as possible. As a further simplification, we assume the volume to be 
cylindrically symmetric around the beam (z)-axis. Thus, the volume is characterized by the 
longitudinal extension L{t) and the transverse radius R{t) and we find 

V{t) = 'rL{t)R\t). (18) 

In order to account for collective flow effects, we boost individual volume elements ac- 
cording to a position-dependent velocity field. For the transverse flow, we make the ansatz 

r]T{r,T)=r/Rr^,{T)r^rf^'{T) (19) 

where Rrmsij) denotes the root mean square radius of the fireball at r and rff^^ij) the 
transverse rapidity at Rrms- 

For the longitudinal dynamics, we start with the experimentally measured width of the 
rapidity interval of observed hadrons 2rij^°"''^ at breakup. From this, we compute the lon- 
gitudinal velocity of the fireball front at kinetic freeze-out vj^°"'^. We do not require the 
initial expansion velocity vl^°^^ to coincide with v(^°"''^ but instead allow for a longitudinally 
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accelerated expansion. This implies that during the evolution 77 = 77^ is not valid (with t]s 
the spacetime rapidity rjg = l/21n((t + z)/(t — z)) unlike in the non-accelerated case. 

The requirement that the acceleration should be a function of r and that the system 
stays spatially homogeneous for all r determines the velocity field uniquely if the motion 
of the front is specified. We solve the resulting system of equations numerically [4^. We 
find that for not too large rapidities 77 < 4 and accelerations volume elements approximately 
fall on curves const. = \Jt'^ — and that the flow pattern can be approximated by a 
linear relationship between rapidity ri and spacetime rapidity r^^ as 77(775) = C^s where C = 
^front I ^fjont ^front -g rapidity of the cylinder front. In this case, the longitudinal 
extension can be found calculating the invariant volume V = J da^u^ as 

(C- 1) ^ ' 

with the spacetime rapidity of the cylinder front. This is an approximate generaliza- 

tion of the boost-invariant relation L{t) = 27/-^''°"*r which can be derived for non-accelerated 
motion. 



5.2. Parameters of the expansion 

In order to proceed, we have to specify the longitudinal acceleration a^(r) (which in turn 
is used to calculate 77{^°"*(r) numerically), the initial front velocity t>Q^°"* and the expansion 
pattern of the radius R{t) in proper time. 

In principle, one would require a = Vp/e. However, our model framework contains a 
homogeneous distribution of matter, therefore Vp = everywhere except at the surface 
of the cylinder. In order to keep this approximation but nevertheless use a more realistic 
acceleration, we make the additional assumption that in a realistic situation a drop in 
temperature would leave the shape of the pressure distribution rather unchanged while 
reducing the overall magnitude. With this assumption, Vp ~ c ■ p. Therefore, we make the 
ansatz 

a. = c. ■ 44 (21) 

e(r) 

which allows a soft point in the EoS where the ratio p/e gets small to influence the acceler- 
ation pattern. Cz and f q^°"* are model parameters governing the longitudinal expansion and 
fit to data. 

Since typically longitudinal expansion is characterized by larger velocities than transverse 
expansion, i.e. ^ w^*^""*, we treat the radial expansion non-relativistically. We assume 
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that the radius of the cyhnder can be written as 

R{t) = Ro + Ct r dr' r dr"^ (22) 
J TO J TO e(r' ) 

The initial radius Rq is taken from overlap calculations. This leaves a parameter ct 
determining the strength of transverse acceleration which is also fit to data. The final 
parameter characterizing the expansion is its endpoint given by tj, the breakup proper time 
of the system. 



5.3. Thermodynamics 

We assume that entropy is conserved throughout the thermalized expansion phase. There- 
fore, we start by fixing the entropy per baryon from the number of produced particles per 
unit rapidity (see e.g. I^l"]). Calculating the number of participant baryons (see [2^) we 
find the total entropy 5*0. The entropy density at a given proper time is then determined by 
s = So/Vir). 

We describe the EoS in the partonic phase by a quasiparticle interpretation of lattice 
data which has been shown to reproduce lattice results both at vanishing baryochemical 
potential /i^ and finite fiB [13] (see these references for details of the model). 

For the phase transition temperature, we choose Tq = 170 MeV based on lattice QCD 
computations at finite temperature for case of two light and one heavy quark flavour (42j |. 
We also note that no large latent heat is observed in the transition and model the actual 
thermodynamics as a crossover rather than a sharp phase transition. Nevertheless, in the 
calculation we assume quarks and gluons as degrees of freedom above Tc and hadrons below 
to simplify computations. Since the time the system spends in the vicinity of the transition 
temperature is small compared with the total time for dilepton emission however, any error 
we make by this assumption is bound to be small as soon as we consider the measured rates 
which represent an integral over the time evolution of the system folded with the emission 
rate. 

Since a computation of thermodynamic properties of a strongly interacting hadron gas 
close to Tc is a very difficult task, we follow a simplified approach in the following: We cal- 
culate thermodynamic properties of the hadron gas at kinetic decoupling where interactions 
cease to be important. Here, we have reason to expect that an ideal gas will be a good 
description and calculate the EoS with the help of an ideal resonance gas model. Using the 
framework of statistical hadronization ^3], we determine the overpopulation of pion phase 
space by pions from decays of heavy resonances created at Tc and include this contribution 
(which gives rise to a pion-chemical potential of order /x^ ^ 120 MeV into the calculation. 
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We then choose a smooth interpolation between decophng temperature Tf and transition 
temperature Tc to the EoS obtained in the quasiparticle description. This is described in 
greater detail in 2^. 

With the help of the EoS and s(r), we are now in a position to compute the parameters 
p{T),e{T),T{T) as well. Since the ratio p{T)/e{T) appear in the expansion parametrization, 
we have to solve the model self-consistently. 



5.4. Solving the model 

In order to adjust the model parameters, we compare with data on transverse momentum 
spectra and HBT correlation measurements. This is discussed in greater detail in 

HQ- 



In j45|, a very similar model is fit to a large set of experimental data, providing differ- 
ent sets of parameters Tf,v±f,Rf,rij^°"'''. Although we use a different (box vs. Gaussian) 
longitudinal distribution of matter, we use the parameters from this analysis as a guideline 
for our transverse dynamics where this difference should not show up and determine r]j^°"'* 



separately. Specifically, we use the set bl from |45[ for the transverse dynamics. 

By requiring R{Tf) = Rj and fy*^""* = f_L/ we can determine the model parameters ct 
and Tf. Cz is fixed by the requirement 7]^'^°'^^{Tf) = rjj^""'^. The remaining parameter Vq^°"'* 
now determines the volume (and hence temperature) at freeze-out and can be adjusted such 
that T{Tf) = Tf. 

The model for 5% central 158 AGeV Pb-Pb collisions at SPS is characterized by the 
following scales: Initial long, expansion velocity fQ^°"* = 0.5c, thermalization time tq = 1 
fm/c, initial temperature Tq = 305 MeV, duration of the QGP phase tqgp = 7 fm/c, 
duration of the hadronic phase Thad = 9 fm/c, total lifetime t/ — tq = 16 fm/c, r.m.s radius 
at freeze-out i^j*"* = 8.55 fm, transverse expansion velocity v±f = 0.537c. 

For the discussion of dileptons, we require the fireball evolution for other than 5% central 
collisions. In this case, we make use of simple scaling arguments based on the initial overlap 
geometry and the number of collision participants. For a detailed description, see [2^0]- 

In j23|, it has been shown that this scenario is able to describe the measured spectrum 
of low mass dileptons, and in Q] it has been demonstrated that under the assumption of 
statistical hadronization at the phase transition temperature Tc, the measured multiplici- 
ties of hadron species can be reproduced. In the model has been shown to describe 
charmonium suppression correctly. None of these quantities is, however, very sensitive to 
the detailed choice of the equilibration time tq. Therefore, we have only considered the 
'canonical' choice tq = 1 fm/c so far. The calculation of photon emission within the present 
framework provides the opportunity to test this assumption and to limit the choice of tq. In 
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|47|, this has been investigated in some detail. Within the present framework, the hmits 0.5 
fm/c < To < 3 fm/c could be found. Variations within these limits, however, do not affect 
the spectrum of dileptons with invariant mass below 1 GeV significantly. 



6. DILEPTON EMISSION 

The emission of dileptons in the model is calculated using the differential emission rate 

lmU{q,T) = -, (23) 



where a = e'^/An and we have neglected the lepton masses. Eq. ()23|) is valid to order a in the 
electromagnetic interaction and to all orders in the strong interaction. Its main ingredient 
is the temperature-dependent spectral function R{q,T). 

The differential rate of eq.()23|) is integrated over the space-time history of the collision 
to compare the calculated dilepton rates with the CERES/NA45 data taken in Pb-Au 
collisions at 158 AGeV (corresponding to a cm. energy of ^/s ~ 17 AGeV) and 40 AGeV 
{^/s ~ 8 AGeV). The CERES experiment is a fixed-target experiment. In the lab frame, 
the CERES detector covers the limited rapidity interval rj = 2.1 — 2.65, i.e. At] = 0.55. We 
integrate the calculated rates over the transverse momentum px and average over rj, given 
that d'^p = MpT dM d-q dpx dO. The formula for the space-time- and p-integrated dilepton 
rate hence becomes 

Tf OO 

d^N 2itM f ^ f rr., f , dN{T{T),M,7],pT) 



dMdr] Ar] 



I dr I dv V{v,TiT)) I dpT PT "^^^^^dtd^p"^^^^^ Acc{M,r^,pr), (24) 

TO 

where Tf is the freeze-out proper time of the collision, V^(?7, T(r)) describes the proper time 
evolution of volume elements moving at different rapidities and the function Acc{M,ri,pT) 
accounts for the experimental acceptance cuts specific to the detector. At the CERES 
experiment, each electron/positron track is required to have a transverse momentum p^ > 
0.2 GeV, to fall into the rapidity interval 2.1 < t] < 2.65 in the lab frame and to have a pair 
opening angle Qee > 35 mrad. Finally, for comparison with the CERES data, the resulting 
rate is divided by dNch/drj, the rapidity density of charged particles. 

As in |23|, we also include the effects of the overpopulation of the pion phase space due 
to decay processes of heavy resonances to the emission from hadronic matter by introducing 
a temperature-dependent pion chemical potential fin(T). 

In addition to the thermal emission of dileptons, we also consider dileptons from vacuum 
decays of vector mesons after the thermal decoupling of the fireball and hard dileptons from 
initial Drell-Yan processes. For details of the calculation of these contributions see ref. 
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FIG. 1: Photon spectral function at T=100 MeV. 



A final remark: The spectral function considered in the present work is only available 
below 1 GeV invariant mass range. Therefore, we do not include thermal emission of dilep- 
tons from hadronic matter above that scale when discussing the Walecka model results. 
This leads to small differences between the two approaches for large invariant masses. This 
region, however, is not dominantly filled by emission from the hadronic phase but by the 
QGP and Drell-Yan contribution. 



7. RESULTS AND DISCUSSIONS 

We now discuss the results obtained for the photon spectral function in the hadronic 
matter arising due to medium modification of the vector mesons with the quantum correction 
effects from baryons and scalar mesons, and its effect on the dilepton emission spectra. 
Contrary to the mean field approximation, there is large drop of the vector meson masses 
due to vacuum polarisations from the nucleon sector. This shifts the p and uo peaks in the 



dilepton spectra to lower values [22, |4^. The additional quantum correction effects from 
the sigma mesons lead to considerable increase of the uj decay width Q]. 

In the present investigation, we choose the renormalised sigma meson self interaction 
coupling \r to be 5, corresponding to the value of incompressibility of nuclear matter to 

are studied in a mixed scenario 
of QGP and hadronic matter using the fireball model as described in the previous section. 
The photon spectral function as obtained in the present hadronic scenario is compared with 
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FIG. 2: Photon spectral function at T=140 MeV and T=160 MeV. 
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FIG. 3: Photon spectral function at T=100 MeV and at densities 0.2po and 0.6p0) with po=0.17 
fm~^. 



the earlier calculation obtained in a chiral SU(3) model. This is described in Il7|l for finite 

n n 

density and 18] for finite temperature modifications. In j22l|, we assumed that it is possible 
to factorize these contributions. 

The spectral functions in both the models are shown in figure Q for vanishing baryon 
density and a temperature of 100 MeV. As expected due to absence of any pronounced 
medium effects for this temperature, the two models yield similar results. 

In figure El the temperature dependence of the photon spectral function is shown. The 
spectral function in the present hadronic scenario is seen to develop a distinct p peak at 
higher temperatures. This is a refiection of the fact that with quantum correction effects, 
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FIG. 4: Dilepton emission rate for 158 AGeV at SPS for the Walecka model with and without the 
colhsion effects. 

the p-mass, due to the tensorial coupUng, has a larger drop as compared to the mass of 
the a;-meson. In the present calculations, the effect of temperature on the p-mass is rather 
moderate up to around a temperature of 160 MeV or so. This is in line with the previous 
calculation of [4^. The broadening of the p peak at higher temperature takes place as the 
Bose enhancement dominates over the effect of drop in p-mass, leading to an increase in 
the p decay width. On the other hand, the density dependence of the p mass is seen to be 
quite significant j22] which leads to the p peak position shifted to lower values in the present 
calculation as illustrated in figure El The mass of the p- meson remains almost unchanged in 
the chiral model, whereas the p- decay width has appreciable enhancement in the medium 
due to inelastic processes, like pN — > ttA^, pN ttA, and also, from pN —>■ uN at higher 
energies. This basically leads to the p peak to be completely dissolved in the chiral model 
jlT^ . leaving only a broad continuum. In the absence of coUisional effects, the decay width 
for p arises from the process p ^ tttt. Here, the Bose- Einstein factors in (0) have the effect 
of increasing the width in the thermal medium, whereas the stronger dropping of the p- 
mass in the medium at higher densities in the present calculations has the opposite effect 
of decreasing the decay width and overcompensates the effect from the Bose enhancement 
factor. The coUisional effects lead to considerable flattening of the spectral function for the 
p peak. 
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FIG. 5: Dilepton emission rate for 40 AGeV at SPS for the Walecka model with and without the 
colHsion effects. 

Figures m and El show the dilepton spectra using the fireball model in the present hadronic 
scenario and a quasiparticle picture for QGP. The inclusion of collisional effects due to 
scattering off by the nucleons is seen to lead to considerable broadening of the p-peak 
in the dilepton spectra. The difference is seen to be more prominent for SPS 40 AGeV 
corresponding to a higher baryon density. The present hadronic scenario leads to a distinct 
structure in the dilepton emission rate due to the undissolved p peak in the absence of the 
contribution from the pN scattering processes, which is seen to be flattened due to these 
collision processes. The dropping of the vector meson masses in QHD in the relativistic 



Hartree approximation |48[ also has been seen to yield similar results. In the absence of 
collisional processes, the drop in the p mass leads to enhancement in the dilepton yield 
below the vacuum p mass, but does not have the rather flat structure as observed in the 
data. Inclusion of inelastic processes pN ttN, pN ttA, pN ujN, give rise to 
considerable broadening of the spectra below the vacuum p-mass. 

Figures ini and m compare the dilepton emission rates of the present investigation to the 
results of chiral model description for hadronic matter as well as to the results from SPS, 
158 AGeV and SPS, 40 AGeV 30% central Pb+Au collision. It may be noted that the fire- 
ball evolution, detector acceptance, contributions from the QGP and Drell-Yan background 
were modelled exactly in the same way for both theoretical curves. The previous hadronic 
scenario based on a chiral model [3| has the observed dissolved p peak arising due to the 
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FIG. 6: Dilepton emission rate for 158 AGeV at SPS for the chiral and Walecka models. 



Bose enhancement factors as well as the inelastic channels. In the present calculations, the 
collisions by nucleons lead to significant broadening of the p-peak in the dilepton spectra. 

Contrary to what has been done in [2^, we have not folded the results with the finite 
energy resolution of the detector in order to indicate what could be observed if the energy 
resolution were increased. The sharp peaks of uj and (p seen in the plots (absent in the 
plots in and seemingly in contradiction with the data) represent vacuum decays of these 
mesons left after kinetic decoupling of the fireball. As apparent from j^l, if the strength 
contained in these peaks is smeared out to represent finite detector resolution, agreement 
with the data is achieved. 

In 0, dilepton and photon emission are studied using a more detailed description of 
the fireball matter using coarse-grained UrQMD and hydrodynamics. The essential scales 
of the expansion dynamics however appear similar to our approach - we also see dilepton 
emission for matter with average temperatures in the region between 120 and 250 MeV for 
a time period of about 15 fm/c. This is no surprise since essential scales of the evolution 
are dictated by the hadronic momentum spectra (and HBT radii) to which our model is 
adjusted. We might therefore expect that simplifying assumptions (such as introducing an 
average T for given r) do not play a large role when one integrates over the complete four 
volume of the expansion. 

It seems that the main difference which might account for the fact that our approach 



18 




Invariant mass [GeV] 



FIG. 7: Dilepton emission rate for 40 AGeV at SPS for the chiral and Walecka models. 



describes the data well between 0.3 and 0.6 GeV invariant mass whereas j50| is somewhat 
on the low side is the use of different spectral functions. However, in order to test this a 
calculation in our framework using yet other spectral functions needs to be carried out which 
is beyond the scope of the present paper. 

In js^], different approaches to describe the dilepton data for 158 AGeV beams are com- 
pared. Unfortunately, each of the calculations shown employs a different fireball model and a 
different description of in-medium modification of the vector mesons, emphasizing the need 
for studies along the line of the present paper. However, the emerging picture for BUU, 
transport and hydrodynamical calculations seems to be that to first approximation, the nor- 
malization of the dilepton invariant mass spectrum is determined by the fireball evolution, 
which is in turn strongly constrained by the hadronic momentum spectra and 2-particle 
correlations and should lead to similar essential scales of the fireball once the builtup of a 
pion chemical potential due to the decay of heavy resonances created at Tc is taken into ac- 
count. The shape of the spectrum is then dictated by the spectral function. Here, a scenario 
employing a dropping in-medium p mass is shown to lead to similar results as broadening of 
the p peak due to the fact that a p mass dropping as a function of T sweeps through a large 
invariant mass region as T changes throughout the fireball evolution. In most approaches 
however, baryons seem to be crucial for the description of the spectral shape, in line with 
our own findings. 
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Unfortunately, the relations of the thermal fireball descriptions investiagted in [51[ to 
essential evolution scales is less clear — for example, no comment is made how transverse 
and longitudinal flow are implemented and what the relation of these models to the actual 
hadronic spectra at freeze-out would be and to what degree the models are tuned to repro- 
duce dilepton data. Nevertheless, the resulting shape of the dilepton spectrum is expected 
to be close to a more detailed calculation whereas the normalization can be expected to be 
less certain. 

Based on the present data, no clear distinction can be made between different approaches 
to calculate the photon spectral function at finite temperature and density. However, choos- 
ing a model which can be shown to describe other hadronic observables as well and using 
different spectral functions within one evolution model seems crucial to make progress in 
understanding the differences between the approaches on a more quantitative level. 



8. SUMMARY 

To summarize, we have investigated here the medium properties of the vector mesons in 
the Walecka model including the vacuum polarization effects and their effects on the dilepton 
spectra using a fireball model. The density dependence of the vector meson properties are 
seen to be the dominant medium effect as compared to the temperature dependence. We 
have considered a mixed scenario of QGP and hadronic matter j23[ and have compared the 
dilepton emission rates of the present calculations to those from a previous calculation using 
a chiral SU(3) model, as well as to the experimental results from SPS at 158 AGeV and 40 
AGeV. 

Since the evolution of the fireball was not in any way adjusted to the dilepton data but 
rather fixed from different observables, we are in a position to compare the different spectral 
functions with the data without being subject to large uncertainties regarding the medium 
evolution. 

Inclusion of vacuum polarisation effects in the Walecka model leads to dropping of these 
masses in the medium, contrary to the mean field approximation. This gives rise to the shift 
of the p and uj peaks in the dilepton spectra to smaller values l^l • Due to the stronger 
medium modification, the p peak is seen to be distinct from the uj peak in the spectra. 
When the scattering due to nucleons is not taken into account in the present calculations, 
there is a more pronounced difference in the two models for SPS, 40 AGeV, as compared to 
SPS, 158 AGeV due to the higher baryon density. However, there is seen to be significant 
broadening of the spectra due to these collision effects within the present investigation. 
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